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Abstract- Based  on  the  classical  Pennes  bioheat  equation  and 
oxygen  diffusion  equation,  a  new  coupled  heat-mass  transfer 
model  for  characterizing  the  cerebral  circulatory  arrest  is 
proposed.  Two  cooling  approaches  (the  surface  cooling  and 
volumetric  cooling)  are  applied  to  analyze  the  effect  of 
hypothermia  on  the  transient  temperature  and  the  oxygen 
consumption  rate  in  different  regions  of  brain.  It  is  shown 
that  the  surface  cooling  contributes  little  to  depress  the 
oxygen  consumption,  which  received  few  attentions  before.  As 
an  alternative,  the  volumetric  cooling  will  significantly  reduce 
the  oxygen  consumption  rate  in  the  brain,  thus  prolong  the 
patients’  survival  time.  The  results,  which  are  close  in 
agreement  with  the  previous  experimental  data,  may  aid  to 
find  more  efficient  treatments  on  the  brain  hypothermia 
resuscitation. 
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I.  INTRODUCTION 

The  brain  is  an  organ  whose  normal  function  depends 
heavily  on  an  uninterrupted  delivery  of  oxygen.  Unlike 
skeletal  muscle  that  can  survive  for  hours  without  oxygen, 
brain  cells  show  irreversible  damage  within  minutes  from 
the  onset  of  oxygen  deficiency.  When  cardiopulmonary 
circulation  is  interrupted  in  many  cardiac  surgical 
procedures  or  accidental  events,  it  will  create  an  imbalance 
between  energy  production  and  consumption  and  cause 
rapid  loss  of  oxygen  and  depletion  of  metabolic  substrate. 
Under  this  adverse  condition,  the  cerebral  protection  by 
hypothermia  has  been  well  established  and  commonly  used 
clinically  due  to  the  decreased  metabolic  requirements  of 
the  cold  brain  tissues.  There  were  extensive  discussion  on 
the  effect  of  hypothermia  and  circulatory  arrest  on  cerebral 
blood  flow  and  metabolism  [l]-[3].  Previous  studies  have 
found  that  total  circulatory  arrest  allows  patients  survive 
for  60  minutes  or  more  at  15-18°C  without  any  subsequent 
cerebral  dysfunction  and  hypothermia  is  effective  for 
protection  against  both  global  and  focal  brain  ischemia. 
Reported  experiments  [4]  suggested  that  the  brain 
continues  to  metabolize  at  a  very  low  basal  rate  during  deep 
hypothermic  circulatory  arrest.  But  to  our  knowledge,  few 
models  have  been  developed  to  quantitatively  investigate 
the  effects  of  cerebral  temperature  changes  during 
circulatory  arrest  to  the  oxygen  consumption  and  how  long 
could  the  patients  withstand  the  oxygen  delivery  deficiency 
until  all  adenosine  triphosphate  (ATP)  stored  are  depleted. 


Therefore,  it  is  necessary  to  develop  a  mathematical  model 
to  address  these  issues  and  explore  a  better  cooling 
approach  to  improve  the  patient’s  ability  to  survive  in  the 
hypoxia. 

In  this  paper,  based  on  the  classical  Pennes  bioheat 
equation  and  oxygen  diffusion  equation,  a  new  coupled 
heat-mass  transfer  model  for  characterizing  the  cerebral 
circulatory  arrest  is  proposed.  Two  cooling  approaches  (the 
surface  and  volumetric  cooling)  are  used  to  analyze  the 
effect  of  hypothermia  to  the  transient  temperature  and  the 
oxygen  consumption  rate  in  different  regions  of  the  brain. 

II .  METHODOLOGY 

The  present  heat-mass  coupling  model  falls  in  the 
appropriate  assumption  that  the  brain  tissue  can  be 
approximately  treated  as  a  uniform  sphere.  The  heat 
transfer  can  be  represented  by  a  one-dimensional  Pennes 
bioheat  equation  in  which  the  local  brain  temperature 
T(r,  t)  is  determined  by 

pC^-  =  K--^^  +  Wb{T,Qm)Cb{E-T)  +  Qm{t,T){\) 
at  r  or 

where  p  ,  C  ,  K  are  density,  specific  heat  and  thermal 
conductivity  of  brain,  respectively;  Cb  is  specific  heat  of 
blood.  These  four  thermal  properties  are  treated  as  uniform 
and  independent  of  temperature.  Qm  is  the  metabolic  heat 
generation  rate;  Wb  is  blood  perfusion  which  may  depends 
on  local  temperature  and  heat  generation  rate;  Ta  is  the 
arterial  temperature. 

At  steady  state  (normal  condition),  the  one  dimensional 
heat  transfer  problem  can  be  established  as: 

K--^P-  +  WbCb(Ta-T)  +  Qm  o=0  (2) 

r  or 

with  boundary  conditions  defined  as: 

-K-^-\r^  =hfo(T-T{0),  r  =  R0  (3) 

—  =  0,r  =  0  (4) 

dr 

where  hm  is  the  heat  convention  coefficient;  T{0  ,  R0 
denote  the  surrounding  air  temperature  and  radius  of  the 
brain  hemisphere,  respectively. 

When  cerebral  circulation  arrests,  a  transient  state,  the 
blood  perfusion  rapidly  reduces  to  zero  in  short  time,  but 
the  brain  metabolism  still  keeps  working  until  it  runs  out  of 
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the  metabolic  substrate.  Therefore,  the  transient  brain 
temperature  T(r,t)  can  be  described  by  the  following  heat 
transfer  equation: 


pcd_L=K±d^Il 

dt  r  dr 2 


+  QmUj) 


(6) 


with  the  same  boundary  conditions  as  given  above  in  the 
steady  state.  The  initial  condition  is  expressed  as: 

T{r,0)  =  T{r),t  =  0  (7) 


Solving  (6)  satisfying  the  boundary  and  initial  conditions 
yields  the  transient  temperature  distribution  in  the  brain 
tissue. 

But  the  toughest  task  is  to  evaluate  the  metabolic  heat 
generation,  which  depends  on  the  local  brain  temperature 
and  time.  To  better  deal  with  this  problem,  we  aimed  at 
developing  a  new  metabolic  heat  generation  model  in  this 
paper.  The  most  widely  used  approach  in  estimating  the 
metabolic  term  has  been  to  set  it  to  be  the  product  of  the 
oxygen  consumption  and  its  caloric  value.  But  this  is  based 
on  the  assumption  that  oxygen,  and  not  glucose,  determines 
heat  generation  in  tissues,  which  is  improper  for  ischemic 
situation.  During  circulatory  arrest,  the  metabolic  rate 
decreases  with  that  of  the  local  temperature.  It  is  assumed 
that  the  metabolic  rate  changes  according  to  the 
temperature  coefficient  (j)  .  On  the  other  hand,  the 
metabolism  dose  not  end  until  the  substrate  was  depleted. 
Then  one  can  assume  that  it  decreases  exponentially  with 
the  time.  Therefore,  the  metabolic  rate  can  be  established  as 
follows: 

Qm (r,  t, T)  =  Qm0  ■  ^nr,o-mo .  exp(_  t/To)  (8) 


where,  Qm0  is  the  reference  metabolic  rate  at  37  °C ,  t0  is 


reference  time  at  37°Cand  is  estimated  as  300  seconds. 
According  to  references  [1],  [2]  ,[5],  [6],  fzl  is  set  as  3.0. 

During  the  circulatory  arrest,  the  transient  oxygen 
behavior  is  governed  by  the  reaction-diffusion  equation: 

(9) 

at  r  or~ 

with  the  boundary  condition  defined  as  dCjdr  =  0  .  r  =  R0 
and  dCjdr  =  0  ,  r  =  0  .The  temperature  dependent  oxygen 
diffusion  coefficient  Dt  is  obtained  from  [7]  as: 

D,(T)  =  a  -exp[Z)-r(r,t)]  (10) 

where  a  =  3.95xl0~10m2 /s  is  a  constant,  b  =  4.6%/°C 
is  the  temperature  coefficient  for  diffusion. 

Commonly  the  oxygen  consumption  can  be  treated  as 
following  the  first-order  kinetics,  i.e.,  M  -  J3C ,  [8]. 
Taking  account  of  the  effect  of  the  temperature,  we  also 
assume  that  oxygen  consumption  varies  with  the  local  brain 
temperature.  Subjected  to  these  assumptions,  the  oxygen 
consumption  rate  is  thereby  expressed  as: 

M(c,,r)  =  ^c,  V^’0-371710  (ii) 

where,  fi  =  4.7  x  1 0  V  is  a  constant,  and  rj  is  temperature 
coefficient  and  set  as  2.5  [1]. 

For  the  calculations,  typical  values  for  the  thermal 
parameters  were  applied  as  given  in  [9]: 

p  =  1.05xl03kg/m3  ,  Cb~C  =  4200J  /  kg  •"  C 


K  =  0.5 W/m  -°  C  ,  Qm0  =  10000W/m3 .  And  R0  =  10cm  is 
used. 

III.  RESULTS  AND  DISCUS  SINS 

Through  numerical  method,  we  calculate  and  analyze  the 
transient  brain  temperature  distribution  and  the  oxygen 
concentration  during  circulatory  arrest.  Under  normal 
condition  (circulatory  arrest  doesn’t  occur),  it  can  be 
regarded  that  the  brain  temperature  is  close  to  37.24  °C  , 
which  is  in  agreement  with  the  physiological  data.  Fig.l  (a) 
shows  the  brain  temperature  distribution  during  circulatory 
arrest  with  heat  convection  coefficient  hf  =  10°C  and 
surrounding  air  temperature  7j  =  20°  C'  .  Clearly,  the 
temperature  near  the  brain  core  (  r  =  0.0  lm  )  increases 
slightly  at  the  beginning,  which  is  mainly  resulted  from  an 
unbalance  between  the  high  metabolic  heat  generation  and 
the  zero  blood  perfusion  cooling.  The  temperature  increase 
in  the  brain  core  at  the  beginning  has  also  been  reported  by 
Xu  et  al.  [5].  The  brain  core  temperature  decreases  much 
more  slowly  than  that  at  the  brain  surface.  Even  after  a  long 
period  of  time  (say  about  120  minutes),  it  still  remains  at 
about  36.15  °C  . 


(a)  surface  air  cooling:  hf  =  10°C  ,  7]  =  20°C 


(b)  surface  ice-water  cooling:  hf  =  30°C  ,  Tf  =  0°C 
Fig.l.  Temperature  distribution  during  circulatory  arrest 
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To  analyze  the  effect  of  hypothermia  on  the  cerebral 
oxygen  consumption  rate,  we  use  two  approaches,  the 
surface  cooling  (water  flowing  around  the  brain  surface) 
and  the  volumetric  cooling  (inserting  cooling  probe  into  the 
deep  brain  or  ventilating  the  cooling  medium  through  the 
mouth  and  nose),  to  examine  the  effect  of  different  cooling 
on  the  temperature  development  and  the  oxygen 
consumption. 


(a)  surface  air  cooling:  h{  =  10°C  ,  T{  =  20°C 


(b)  surface  ice-water  cooling  hf  =  30°C  ,  Tt  =  0°C 
Fig.  2.  The  oxygen  concentration  distribution  during 
circulatory  arrest 


Fig.l  (b)  depicts  the  brain  temperature  distribution 
under  convective  cooling,  h{  =30°C  and  T{  =  0°C  .  The 
temperature  near  the  brain  core  (  r  =  0.01m )  still  increases 
slightly  from  37.24  °C  to  the  maximum  37.73  °C  in  about 
26  minutes,  then  decreases  slowly.  The  temperature  on  the 
brain  surface  can  be  significantly  reduced  by  the 
convective  ice-water  cooling,  but  the  effect  of  this  cooling 
on  the  deep  brain  is  not  evident  since  it  takes  82  minutes  for 
the  brain  core  temperature  to  change  from  37.29  °C  to 
36.5  °C  .  These  results  are  consistent  with  the  previous 
clinical  data.  Surface  cooling  of  head  and  irrigation  of  the 
nasopharynx  with  ice  water  lowered  deep  brain 
temperature  from  37  °C  to34°C  only  after  30-60  minutes 
[5].  Mellergard’s  [10]  and  Olsen  et  al’s  [11]  studies  also 
demonstrated  that  isolated  head  cooling,  whether  with 
frozen  liquid  or  a  cooling  helmet,  had  a  very  limited  effect 
on  lowering  the  brain  temperature.  As  shown  in  Fig. 2,  the 


oxygen  consumption  rate  differs  little  between  the  two 
cases  of  ht  =  10°C  ,  Tt  =  20°C  and  ht  =  30°C  ,  Tt  =  0°C  . 
Overall,  for  surface  cooling,  the  stored  oxygen  amount 
decreases  sharply.  It  is  almost  completely  depleted  within 
less  than  20  minutes.  At  this  time,  the  oxygen  concentration 
is  three  orders  of  magnitude  lower  than  the  initial  oxygen 
concentration.  From  (10)  and  (11),  it  is  seen  that  if  the 
temperature  of  the  brain  can  not  be  quickly  reduced,  both 
the  oxygen  diffusion  coefficient  D,(T )  and  the  oxygen 
consumption  rate  M(C \  T )  will  keep  at  a  higher  value  due 
to  the  weak  cooling  effect.  Furthermore,  the  stored  oxygen 
amount  is  depleted  so  quickly  that  it  can  not  meet  the  basal 
metabolic  requirement,  which  may  cause  the  brain  injury. 
Therefore,  the  surface  cooling  is  in  fact  not  very  useful  in 
reducing  the  oxygen  consumption  and  thus  preventing 
patients  from  the  brain  injure. 


0.15 


0.00 


(a) 


Fig.  3.  The  oxygen  concentration  distribution  during  circulatory 
arrest  using  volumetric  cooling  and  the  uniform  brain 
temperatures  are  assumed  as:  (a)r  =  5°C  (b)7’  =  15°C 

In  order  to  improve  the  brain  hypothermia  resuscitation 
and  prolong  the  patient’s  survival  time,  an  alternative 
cooling  approach  -  the  volumetric  cooling  is  explored.  In 
this  way,  a  significant  temperature  decrease  in  the  brain 
core  can  be  realized  by  ventilating  cold  water  to  patient’s 
circulatory  system  or  directly  inserting  the  cooling  probe 
into  the  deep  brain.  The  oxygen  concentration  distribution 
is  depicted  in  Fig. 3  on  assuming  that  all  the  brain 
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temperature  positions  are  fixed  at  a  uniformly  lowered 
temperature.  Compared  with  the  case  of  surface  cooling, 
the  oxygen  concentration  decreases  much  more  slowly. 
When  the  uniform  brain  temperature  alters  from  15  °C  to 
5  °C  ,  the  average  oxygen  concentration  decreases  by  46% 
and  15%  after  about  20  minutes,  respectively.  Fig. 4  clearly 
shows  that  at  the  region  near  the  brain  core  (r  =  0.01m) 
and  under  uniform  volumetric  brain  temperature  T  =  5"C  , 
the  oxygen  concentration  is  still  at  about  0.0021  kg/m3  at  20 
minutes.  Further  calculations  demonstrate  that  the  oxygen 
consumption  rate  at  this  case  is  several  orders  of  magnitude 
lower  than  that  at  the  case  of  surface  cooling.  Therefore,  it 
is  likely  that  the  temperature  decrease  in  the  whole  brain 
ensures  the  oxygen  demand  sufficient  to  permit  tolerance 
for  a  long  period  of  absent  oxygen  delivery.  In  recent 
studies,  Xu  et  al  [6]  pointed  out  that  the  circulatory  cooling 
may  be  essential  for  a  rapid  brain  cooling  required  for 
cerebral  protection.  Mariak  et  al’s  [12]  experimental  data 
demonstrated  that  heat  loss  from  the  upper  respiratory  tract 
can  directly  affect  the  intracranial  temperature  in  human. 
But  this  method  is  not  practically  used,  which  is  the  issue  to 
be  solved. 


Fig.4.  The  effect  of  two  cooling  approaches  on  oxygen 
concentration  in  the  inner  brain  region  (  r  =  0.01m  ) 

There  are  several  assumptions  and  limitations  in  this 
study  must  be  pointed  out.  First,  it  is  assumed  that  blood 
perfusion  immediately  reduces  to  zero  following  the  onset 
of  circulatory  arrests,  which  is  the  extreme  case  in  the 
present  simulation.  However,  calculations  show  that  there 
is  a  negligible  oxygen  concentration  difference  between  the 
case  of  zero  perfusion  and  that  of  keeping  perfusion 
unchanged  for  short  time.  In  this  model,  it  is  simply 
assumed  that  the  thermal  properties  are  uniform, 
independent  of  space  and  temperature.  In  fact,  the  heat 
generation  is  a  space  dependent  value.  A  spherical 
geometry  is  the  simplest  one  that  describes  the 
one-dimension  temperature  distribution  from  interior  of  the 
brain  to  the  skull.  Therefore,  we  and  other  researchers  [5] 
use  this  symmetrical  geometry  model.  During  circulatory 
arrest,  the  oxygen  consumption  rate  is  proposed  to  take  the 
new  form  as  (11).  According  to  previous  experimental 
study  [1],  the  temperature  coefficient  considerably  varies 
with  the  temperature.  It  falls  in  2. 0-3.0  at  temperatures 
between  37  °C  -28  °C  .  However,  when  temperature  is  less 


than  28  °C  ,  the  reported  value  is  approximately  two  times 
of  that  at  the  higher  temperature.  In  this  paper,  the 
temperature  coefficient  was  treated  as  constant  and  equal  to 
3.0.  However,  although  the  oxygen  consumption  rate  was 
perhaps  overestimated  at  lower  temperature,  the  above 
consideration  does  not  change  the  conclusions 
substantially. 

IV.  CONCLUSION 

Thermal  factor  is  one  of  most  important  factors 
responsible  for  the  oxygen  delivery  deficit.  Results  in  this 
paper  demonstrated  that  the  surface  cooling  may  provide 
some  advantage  by  cooling  the  brain  cortex  peripherally. 
But  its  effect  to  depress  the  metabolic  oxygen  is  very 
limited.  In  contrast,  the  volumetric  cooling  has  the  potential 
to  substantially  depress  the  oxygen  metabolism,  which 
points  out  the  new  way  to  develop  high  efficient  cooling 
system.  Further  experiments  and  clinical  practices  will  be 
performed  in  the  near  future  to  test  our  theoretical 
predictions. 
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